clear all;
x0=[0,0]';
[x,fmin,k] = punishment(x0);%惩罚函数法